

function infected=find_SImodel(adj,i,steps,p)

inf_list=zeros(size(adj,1),1);
inf_list(i)=1;
for ii=1:steps
    candidates=inf_list'*adj;
    for jj=1:size(adj,1)
        if candidates(jj)>0 
            for kk=1:candidates(jj)
                if rand<p
                    inf_list(jj)=1;
                end
            end
        end
    end
end

infected=sum(inf_list);